SHEAR-FREE PERFECT FLUIDS WITH LINEAR 
EQUATION OF STATE 



RADIX SLOBODEANU 



Abstract. We prove that shear-free perfect fluid solutions of Ein- 
stein's field equations must be either expansion-free or non-rotating 
(as conjectured by Ellis and Treciokas) for all linear equations of 
state p = wp except for w E {-f , -\, ±, \). 
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For the Einstein field equations with perfect fluid sources, an impor- 
tant insight into the structure of solutions can be obtained under two 
fairly general and observationally defensible assumptions: the equa- 
q-i tion of state for the pressure p and energy density p is barotropic (i.e. 

p = p(p)) and the flow is shear-free (i.e. fluid's velocity is a transver- 
ki); sally conformal vector). While the former restriction is specific to isoen- 

tropic fluids, the latter is essentially a kinematic condition of isotropy. 
In cosmology, the shear-free condition expresses the isotropy of the 
relative recessional motion of the galaxies (but allows the red shift 
O | and the cosmic microwave background radiation to be anisotropic) and 

is a common feature of standard spacetimes, such as the FRW and 
Godel models. In kinetic theory this characterizes [35J the velocity of 
q ! collision-dominated gases with isotropic distribution function, under 

CO ; the Einstein-Boltzmann equations. 

The first result indicating the nature of solutions in this context was 
stated by Godel (|12j): a shear- free dust fluid (i.e. p = 0, in particular a 
^ ! congruence of timelike geodesies) in a spatially homogeneous spacetime 

of type IX cannot both expand and rotate. This contrasts Newtonian 
cosmology where expanding and rotating (shear-free) solutions exist 
and can avoid the singularity formation [21]. Godel's result was later 
proved [T5] to remain true without symmetry assumptions and it can be 
seen as a timelike analogue of the Goldberg- Sachs theorem [TS]. After 
proving that the same conclusion holds for radiation fluids (i.e. p = 
Ellis and Treciokas Q35J) formulated the shear-free fluid conjecture: 
If the velocity vector field of a self-gravitating barotropic 
perfect fluid (p + p ^ and p = p(p)) is shear-free, then 
either the expansion or the rotation of the fluid vanishes. 
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Since then, many partial results have been obtained by specializing 
(i) the curvature (divergence-free electric or magnetic part of Weyl cur- 
vature [S1EZ], Petrov types N [3] and III [I]), (ii) kinematical quantities 
(parallel vorticity and acceleration [38J , functionally dependent expan- 
sion and energy density [20J or expansion and rotation scalar [33J , fluid 
flow parallel with a conformal vector field j8]), or (Hi) the equation of 
state (p = wp, with w = 0, ±~, § [13 HH1 [321 E5J EH]). Further aspects 
supporting the conjecture can be found in [10J; for a recent discussion, 
see [T6j . 

In Riemannian geometry, several results of similar flavour have been 
independently discovered. They are all related to the notion of (r-) 
harmonic morphism (see pQ for a general account). To connect them 
to the above results in general relativity, let us recall that the relativis- 
tic perfect fluids (whether self-gravitating or not) can benefit from a 
variational treatment if we define the action functional in terms of a 
submersion from the (4-dimensional) spacetime having its fibres tan- 
gent to the fluid velocity, U . This approach, initiated in [6], [TTJ [33] , has 
been mathematically settled in [7] and renewed in recent physics litera- 
ture, e.g. [131 HI]. in particular, it turns out [31] that a shear-free rela- 
tivistic perfect fluid with linear equation of state (U,p = r j^-p,p = A r ) 
(r 7^ 0) is determined by a (locally defined) horizontally conformal 
submersion which is a critical point of the action functional, i.e. an r- 
harmonic morphism (|21[) with 1-dimensional timelike fibres tangent to 
U and dilation A. In Riemannian signature and dimension at least four, 
precisely in the 1-dimensional fibres case, 2-harmonic morphisms on 
Einstein manifolds and (r > 2)-harmonic morphisms on constant cur- 
vature spaces have been classified ([2JI22JI2S1I2Z]): they are all either of 
Killing type (with U collinear with a Killing vector field and div U = 0, 
so " expansion- free" ) , of warped product type (with U self-parallel and 
orthogonal to a foliation by hypersurfaces, so "non-rotating"), or (only 
on Einstein four-manifolds) of a third type which implies Ricci flatness. 
A similar result holds also on conformally-flat manifolds [28J. 

Adopting the standpoint of r-harmonic morphisms, in this paper we 
prove the shear-free fluid conjecture for all linear equations of state 
p = wp except for six values of w which require a case by case analysis. 
Together with the four values already settled, they form the full set of 
exceptional values of the present proof. 

2. Preliminaries 

In this section we present the terminology related to perfect fluids 
and introduce some basic results needed latter on. For a summary of 
notations and conventions, see the Appendix. 

2.1. Perfect fluids coupled with gravity. 
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Definition 1. ([23, 25J) Let (M,g) be a 4- dimensional spacetime. A 
triple (U,p,p) is called (relativistic) perfect fluid if: 

(i) U is a timelike future-pointing unit vector field on M , called the 

flow vector field (or normalized 4-velocityj, 
(ii) p, p : M — > R are real functions, called mass (energy) density 

and pressure, respectively, 
(in) the stress-energy tensor of the fluid is conserved: 

div (pg + (p + p)ru ® w) — 0, 

where w(X) = g(U,X), for all X tangent to M. 
If instead of (Hi), the Einstein field equations are satisfied 

(1) Ric — | Seal -g = p g + (p + p)vj ® w, 

then (U, p, p) is called perfect fluid coupled to gravity ( or self-gravitating/ 

Condition (Hi) of the definition splits into the fluid's equations: 

(p + p)VjjU + grad^p = (Euler equations), 
(p + p) div U + U(p) = (energy conservation), 

where grad w p is the spatial pressure gradient. The coupling with grav- 
ity imposes the following block- diagonal form of the Ricci tensor: 

Ric(U, U) = \(p + 3p), Ric(X, U) = 0, VX _L U, 

(3) Ric(X,Y) = l(p-p)g(X,Y), VX,Y±U. 

Notice that Seal = p — 3p. 

We denote by V the (vertical) foliation spanned by U and by % its 
(horizontal) complement in the tangent bundle of M. The orthogonal 
projection on TL (or V) of a vector field X will be denoted by X n (or 
X v ). Then any vector X has the splitting: X = X n — zu(X)U. Let 
W be the tensor defined by (VU)(X, Y) = g(V x U, Y), for all X, Y. 

Proposition 1. The covariant derivative of the flow vector field de- 
composes as follows 

(4) VU + w® (V v U) b =u + o- + ±divUg H , 

where g n (X,Y) = g(X n ,Y n ) is the orthogonal projector to U and 
cu(X,Y) = ±dw(X n ,Y n ), 
a(X, Y) = \ (C u9 - § div U g) (X n , Y n ). 

Proof. Notice that Koszul formula can be rewritten as follows: 

g(V x U, Y) = dw(X, Y) + \(C u9 )(X, Y), VX, Y. 

Since g(V x U, Y) = g(V x uU, Y n ) - w(X)g(V v U, Y), applying Koszul 
formula on the first right-hand term yields the result. □ 
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Definition 2 (Kinematical quantities). For a relativistic fluid (U,p, p), 
Vfy£7 is the acceleration (vector), divU is the expansion rate (scalar), 
bj is the rotation rate or vorticity (2- form) and a is the shear rate 
(symmetric and traceless 2-covariant tensor). 

From the following slightly reformulated Yano's identity (|40j) 

(6) Ric(X,Y) =div(V x Y)-X(dwY)-l(C x g,C Y g) + ^(dX\dY b ), 

where X, Y are tangent vectors to a (semi-)Riemannian manifold, we 
obtain 

Proposition 2 (Raychaudhuri equation [25]). On the spacetime (M, g) 
the following identity hold: 

(7) U(divU) = -§(div[/) 2 + 2 (M 2 - H 2 ) +div(V c/ f/)-Ric(f/,[/). 

In this context, the shear-free fluid conjecture reads: if Equation 
(jTJ is satisfied by a fluid with p = p(p), p + p ^ and a = 0, then 
(div U)u = 0. We shall treat the special case described below. 

2.2. Shear-free linear-barotropic fluids. Now consider a shear-free 
(a = 0) perfect fluid with linear equation of state p = wp (which is 
the 'y-law used in cosmology). We can always write it in the form 
(U,p = ^A r ,p = A r ), where A is some positive function (locally) on 
M (A -1 represents the length scale) and r e E*. 
The fluid equations (j2J) become, respectively 

(8) - Vt/£/ + (3-r)grad w (lnA) = 0, 

(9) div U + 3U(\n A) = 0. 

Remark 1. According to [2j [22], (J8]) and (J9j) are the necessary and 
sufficient conditions for U to be tangent to the fibres of a (locally 
defined) r-harmonic morphism ip : (M,g) — > (N,h) into some 3- 
dimensional Riemannian manifold, with dilation A. Recall that [2T] an 
r-harmonic morphism ip is characterized as a horizontally conformal 
(i.e. ip*h = X 2 g n ) critical point of the action functional j M |dv?| r vol 9 . 
The corresponding Euler- Lagrange equations are given precisely by (jSJ), 
while d2|) is an identity fulfilled by any horizontally conformal submer- 
sion (for more details see [3T] ). 

In this case, the Conditions <^ for the fluid to be coupled with 
gravity turn out to be 

Ric(U, U) = ^A r , Ric(X, U) = 0, VX _L U, 

^ Ric(X, Y) = ^f\ r g(X, Y), VX, Y _L U. 

and therefore Seal = (4 — r)A r . 

A crucial notion in the sequel are basic tensors. 

Definition 3. ([lj) A section q of ((g> r "H) ® (® S "H*) is called basic if 
(C w q) H = for allW G T(V). 
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In particular, a function / on M is basic if W(f) = 0, and a hori- 
zontal vector field X on M is basic if [W,X} H = 0, for all W G T(V); 
the latter condition means that X is projectable on N. 

Let us define the fundamental vector field of V by 

V = X 3 ~ r U. 

Analogously to [26] (see also [H p.341]), we have 
Proposition 3. The Euler equation (jSj) is equivalent to 
(11) [V,X} = 0, 

for all basic vector fields X on M. 

Proof. Let X be a basic vector field. Then [V, X] is vertical. Since 

g([V,X],U) = \ 3 - r g{(3 - r)grad w lnA - V V U,X), 

we conclude that [V, X] = if and only if Euler equation is satisfied. □ 

Let be the 1-form dual to V, defined by d(X) = -X 2 ^ g(X, V) 
for all X, and let Q = d-d. 

Lemma 1. Let (U,p = £^X r ,p = X r ) be a shear-free perfect fluid. 
Then 

(i) Sl(X, Y) = \ 2r - 6 g([X, Y],V) for all X, Y e T(H); 

(ii) Q = —2\ r ~ 3 ui (Q is proportional to the rotation ratej; 
(Hi) = 0, for all W G T(V) (Q is a horizontal 1-form) 

(iv) Q = if and only if % is integrable Q; 

(v) £w^ = or, equivalently, W(Q(X,Y)) = 0, for all X,Y basic 
vectors and for all W G T(V) (Q is a basic 1-form). 

(vi) For any basic vector field X, the following function is basic: 

\- 2 [5Q(X) + (r - 4)fi(X, grad In A)] . 

Proof, (i) and (ii) are immediate using only the definitions. For (Hi) 
we can check that Q(V,X) = for a basic vector X by employing 
Equation (TTTT) . (iv) is the Frobenius' theorem upon applying (Hi), (v) 
is the consequence of Cartan's formula Cw^ = iwdQ + d(iwty- 

To prove the last assertion (vi), let {X,Y,Z} be orthogonal basic 
vectors from a preferred frame (see Remark below). By direct compu- 
tation we obtain 

5Q(X) = A 2 [(4 - r)Y (In X)Q(X, Y) + Y(Q(X, Y)) + X 2 g([Z, Y], Z)Q(X, Y)] , 
5Q(Y) = X 2 [(4-r)X(ln\)tt(Y,X) -X(Q(X,Y)) - X 2 g([Z,X], Z)Q(X,Y)] , 
5Q(Z) = -A 4 c/([X,F],Z)fi(X,F). 

Since [X, Y] H , [X, Z] and [Y, Z] are basic vectors, the functions X 2 g([X, Y],Z), 
X 2 g([Z, X], Z) and X 2 g([Z, Y], Z) are basic. Using also the fact that Q 

1 In other words, for the 4-velocity of a perfect fluid, being hypersurface orthog- 
onal is equivalent with being vorticity-free (i.e. having zero rotation). 
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is basic cf. (v), we obtain the conclusion for X, Y and Z, and therefore 
for any basic vector. □ 

Remark 2 (Preferred orthonormal frames, cf. [H|26]). We can always 
choose a local orthogonal frame {X, Y, Z} of basic horizontal vector 
fields and we may suppose that their lengths satisfy \X\ = \Y\ = \Z\ = 
1/A and that Z satisfies tz^i = 00 (i.e. Z is collinear to the vorticity 
vector {*^uj)^). Then the contractions andiy^l are both basic and 
orthogonal, i.e. (ixQ,iyQ) = 0. Since dQ(X,Y, Z) = 0, we also have 

(12) (C z n)(X,Y) = 0. 

Moreover, \Q\ 2 = X A Q{X,Y) 2 and \i x tt\ 2 = \tytt\ 2 = X 2 tt(X,Y) 2 . 
We call {U, XX, AY, XZ} a preferred orthonormal frame □. 



3. Constraint equations 

The very existence of a shear-free perfect fluid constrains the space- 
time geometry, in particular the Ricci curvature (Proposition H]) that, 
in addition, must have block- diagonal form ([3]) by the Einstein field 
equations. For a shear-free fluid with linear equation of state p = 
with p = X r and r ^ 2, we see these curvature restrictions as constraints 
on the second order derivatives of In A. This will provide compatibility 
conditions at the level of 3 rd order derivatives, in the form of polynomial 
relations in first derivatives only. 

In this section we collect the constraint equations that are useful in 
the proof of the conjecture. Analogously to (see also [U p. 343]) we 
have 

Proposition 4 (Ricci curvature restrictions). Let (M,g) be a four di- 
mensional spacetime and if : (M, g) — > (N, h) an r-harmonic morphism 
of dilation X, into a 3- dimensional Riemannian manifold. If the fibres 
are tangent to the timelike unit vector field U , then the following iden- 
tities hold for all horizontal vectors X, Y : 
(13) 

Ric(U,U) = (3-r)AlnA+(6-r)f/(t/(lnA))+3(r-4)f/(lnA) 2 +^^|fi| 
(14) 

Ric(X, U) = 2X{UQn\)) - ^{5Q{X) + 2(3 - r)Q{X, gradln A)}, 

Ric(X, Y) = <p* Ric N (X, Y) + (r - 2) Hessi nA (X, Y) 

(15) - [(r - 3) 2 + l]X(ln X)Y (In A) + ^f^{i x ^ ly®) 

+ (A In A + (r - 2)| gradln A | 2 )g(X, F). 



Since Q, is basic, it locally descends to a 2-form on N that can be seen as a skew- 
symmetric linear map TN — > TN. Since dim is odd, such a map is singular. 
3 A similar choice of an orthonormal tetrad was done in [36] . 
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Proof. The first identity is simply another way of writing Raychaudhuri 
equation (J7J) in the shear- free case, by using Equations (jSJ) and (J5J). 

The second identity can be derived directly from Yano's formula ([H]), 
taking into account the decomposition (jJJ). This formula is known in 
physics literature as shear-divergence identity or (0, a)— equation and 
it is one of the standard constraint equations, cf. e.g. [30] . 

The third identity is obtained by taking the trace of the following 
formula, true for any horizontally conformal map ip : (M, g) — >■ (N, h) 
and for all horizontal vectors X, Y, Z,T ([1, p. 320]) 

(R(X,Y)Z,T) = X- 2 (R N (dp(X),d V (Y))dp(Z) 1 dp(T))} 

- (X(ln X)Y - Y (In X)X, T(ln \)Z - Z(ln X)T) + {{Y, Z) Hessi nA (X, T) 

- (X, Z) Hess lnA (Y, T) + (X, T) Hess lnA (Y, Z) - (Y, T) Hess ln A (X, Z)} 

+ i{(/(x, z), /(y, r)) - (j(y, z), /(x, t)) + 2(/(x, y), i(z, t))} 

+ ( (Y, Z) (X, T) - (X, Z) ( Y, T) ) I grad In A | 2 , 

where I(X,Y) = [X, Y] v , combined with the propagation of shear: 

R{U, X, 17, Y) + |(£ V ^)(X, Y) - § [C/(div [/) + |(div f/) 2 ] <?(X, Y) 
+ (i x u, iyuj) + ff(Vt,tf, X)g(VuU, Y) = 0, VX, Y e H. 

□ 

Using the gravity coupling condition ffTOl and Remark [Tj we obtain 

Corollary 1. 7ei (Z7,p = ^^X r ,p = A r ) be a shear-free perfect fluid 
coupled with gravity on (M,g). Then the identities ([TBI . (TH1) and (|T5|) 
/io/d u>zi/i i/ie left-hand side replaced by z -^A r , and ^pA r g(X, Y) ; re- 
spectively, where N is endowed with a metric h such that the projection 
along U, p : (M,g) -> (X, h), satisfies <p*h = X 2 g n . 

Taking the (horizontal) trace of f fT5|) and combining with f|T3|) gives: 

Corollary 2 (Trace constraints). On a spacetime (M,g) with a shear- 
free perfect fluid (U,p = r ^X r ,p = A r ) coupled with gravity, the fol- 
lowing identities must hold 

A In A =^U(\nX) 2 - ^X r + ^A 6 - 2r |fi| 2 
(16) " ' 



fr-f- 6 ) | grad^ In A| 2 + ^A 2 Seal", 



U(U(\nX)) = - ^U(\nX) 2 + ^X r - ^A 6 - 2r |fi| 2 

' 1 ' ' _ (2r-5)(r-3) | ^ A |2 _ ^2 g^jTV 



3.1. Constraints in a preferred frame. Let {AX, AY, AZ, U} be a 

preferred orthonormal frame. In the following we will obtain some 
constraints for the existence of a shear-free fluid (U,p = r -^X r , p = X r ) 
coupled with gravity in terms of these basic vector fields X, Y and Z. 
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Condition (TTTj) is equivalent to 



(18) 

X(V(ln A)) =^±{f3{X) + (10 - 3r)Y (In X)Q(X, Y)} + (3 - r)X(ln A) F (In A), 
K(V(ln A)) =^pl{/3(y) _ (10 - 3r)X(ln A)Q(X, Y)} + (3 — r)Y (In A)V(ln A), 
Z(V(lnA)) = ^>/3(Z) + (3 - r)Z(lnA)V(lnA), 

where (3(T) := \~ 2 [5Q(T) + (r — 4)fi(T, grad In A)] is basic whenever T 
is a basic vector, according to Lemma [TJ 

Assume r ^ 2 and let A(r) = {r ~% +1 and B(r) = 1 + {2r ~^ 6) ■ 
Condition ()15j) combined with the trace constraints gives: 

(19) 

Hessi nA (X,X) =(A(r) - 5(r))X(lnA) 2 - 5(r)[F(lnA) 2 + Z(lnA) 2 ] + ±\- 2 U(\n\) 2 

- ^ ( V * Ric N (X, X) + ^ Scal^) , 
the analogous equation for Hessi nA (Y, Y) and 
(20) 

Hess lnA (Z, Z) =(A(r) - B(r))Z(\n\) 2 - B{r)[X (In A) 2 + F(lnA) 2 ] + |A- 2 [/(lnA) 2 

_1_ 1 \r-2 3r-14 \8-2rp/ y y^2 
+ 6 A - 24(^2) A i2(A,Y) 

-^- 2 (ip*Ric N (Z,Z) + r -f£Sc8l N ). 
Condition (TT3T) on pairs of orthogonal vectors gives: 
(21) Hess ln X (X, Y) - A(r)X(\n X)Y (In A) + Ric^X, F) = 

and similar equations for (X, Z) and (Y, Z). 

The previous equations prescribe second order derivatives of In A in 
terms of its first derivatives. Differentiating them along V (i.e. "propa- 
gating" ) and using commutation ( ITT]) will provide us with compatibility 
conditions purely in terms of first derivatives. 

By taking the derivative along V of the equation (cf. ()15p ) 

Hess ln X (X, X) - A(r)X(hx A) 2 + ^ Ric^X, X) 
' 11 ' = Hess lnA (y, Y) - A{r)Y{\n A) 2 + ^* Ric^Y, Y), 
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inserting X(V(lnA)), Y(V(lnX)) from (USD and X(Y(\n\)) from (j2T]) . 
and simplifying the result using again ( I22p . we obtain: 

(3 - r)V{ln A) { (5 - r - A(r)) [X(ln A) 2 - F (In A) 2 ] 
+ ^p* Ric JV (F, Y) - ^* Ric^X, X)] } 
= - ^{2(10 - 3r)(13 - 3r - A(r))ft(X, F)X(ln A)F(ln A) 
+ [(15 - 3r - 2A(r))/3(X) + (10 - 3r)F (ft(X, K))]X(ln A) 
(23) + [-(15 - 3r - 2A(r))/3(Y) + (10 - 3r)X(Q(X, F))]K(ln A) 
+ (10 - 3r)A 2 #(Vx^ + VyX, Z)fi(X, Y)Z(\n A) 
+ X(/3(X)) - Y{(5{Y)) - 2 -^ 1 n{X : Y)ip* Ric N (X, Y) 
+ X 2 g([X, Y],Y)f3(X) + A 2 <?([X, Y],X){3(Y) 
+ A 2 (g([X,Z],X)-g([Y,Z],Y))/3(Z)}. 

Analogously, by propagating the following equation (cf. (TT5]) ) 
(24) 

Hess lnA (X, X) - A(r)XQn A) 2 + ^<p* Ric 7V (X, X) + ^X s ~ 2r Q(X, Yf 
= Hess lnA (Z, Z) - A{r)Z{\n A) 2 + ^ Ric 7V (Z, Z), 

we obtain 
(25) 

(3 - r)V{la A){(5 - r - A(r))[X(ln A) 2 - Z(lnA) 2 ] 
+ ^* Ric^Z, Z) - Ric N (X, X)]} 

= - ^{(10 - 3r)(13 - 3r - A(r))fi(X, F)X(ln A)F(ln A) 

+ [(15 - 3r - 2A(r))0(X) + (10 - 3r)A 2 #([Z, Y], Z)Sl(X, Y)]X(\n A) 

+ (10 -3r)[X(Q(X,Y)) - X 2 g([Z,X},Z)n(X, Y)]YQn\) 

+ [-(15 - 3r - 2A(r))/3(Z) + (10 - 3r)X 2 g(V x Y, Z)£l(X, Y))]Z(\n A) 

+ X(/3(X)) - Z(/3(Z)) - ^n(X, Y)(p* R\ C N (X, F) 

- A 2 s([Z, X], Z)/3(X) + X 2 (g([X, Y], X) - g([Z, Y], Z))f3(Y) 

+ X 2 g([X,Z],X)P(Z)}. 

Notice that, inside the brackets of both (|2"3"|) and (1251) . the polynomial 
expressions in the derivatives of In A have basic coefficients. 

Now we exploit the commutation of the covariant second order hor- 
izontal derivatives. Since [X, Z] is basic, insert ffTSl) in 

X(Z(V(lnX))) - Z(X(y(lnA))) - [X, Z](V(lnX)) = 0, 
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then substitute the 2 nd order derivatives of In A by means of 0211) . to 
obtain 

(r-4)(10-3r) ^ ( ^ Y )Y(\n A)Z(ln A) 

+ {(5 - r)P(Z) - (10 - 3r)fi(V x Z,X)}X(lnA) 

(26) -(10-3r)Q(X,[Z,Y])Y(\nX) 

+ {-(5 - r)/3(X) - (10 - 3r)X 2 g([Z, Y], Z)Q(X, Y)}Z(\n A) 

+ d/3(X, Z) + *%Efn(X, Y)ip* Ric(F, Z) = 0, 

whose left-hand term is a polynomial expression in X(lnA), Y(ln\) 
and Z(\nX) with basic coefficients, denoted for simplicity as follows: 
fey (In A)Z(ln A) + &iX(ln A) + b 2 Y (In A) + b 3 Z(\n A) + b . 

Similarly, from Y(Z(V(\n A))) - Z(Y(V(\n A))) = [Y, Z](V(\n A)) (or 
simply by permuting X and Y in the above relation) we obtain 

- (r ~ 4 ^~ 3r) Q(X,F)X(lnA)Z(lnA) 
+ (10 - 3r)Q([Z, X], F)X(ln A) 

(27) + {(5 - r)0(Z) - (10 - 3r)fi(VyZ, Y)}Y(\n A) 

+ {-(5 - r)/3(Y) + (10 - 3r)X 2 g([Z, X], Z)Q(X, Y)}ZQn A) 

+ d/3(Y, Z) - 1J ^n(X, Y)<f* Ric N (X, Z) = 0, 

whose left-hand term has basic coefficients too and is written down as 
-&X(ln X)Z(\n X) + Cl X(ln A) + c 2 Y (In A) + c 3 Z(ln A) + c . 

Finally, using the fact that [X, Y] n is basic, insert ffT8|) in 

X(Y(V(lnA))) - Y(X(V(lnX))) - [X, Y](V(ln A)) = 0, 

then replace the 2 nd order derivatives of In A by means of fTT7|) and ffT9l . 
to obtain 

n(X,Y){C(r)[X(\nX) 2 + F(lnA) 2 ] + D(r)Z(\nX) 2 } 

(28) + (2iX(ln A) + d 2 Y (In A) + rf 3 Z(ln A) + L(A r_2 , A 8 " 2r ) 
+ 5(4 - r)n(X, Y)X- 2 U(\n A) 2 = 0, 

Where C(r) = _ Hilzfl|i*30) ^ = _ (r-4)(2r-_5)^-18) > 

L(A r_2 , A 8_2r ) = d/3 H (x, y){^A r - 2 - 15r 1 2 - r 5 ! r + 40 A 8 - 2 ^(x, r) 2 + 

i^T [v* ^(X, X) + ^* Ric JV (F, Y)} - (r - 6 y 2 7 18) Seal* }, 
d x = (5 - r)/3(Y) - (10 - 3r)X(fi(X, F)), d 2 = -(5 - r)/3(X) - (10 - 
3r)Y(Q(X, Y)) and d 3 = -(10 - 3r)Z (fl(X, Y)) (notice that the func- 
tions di's are basic). 

4. Strategy of the proof 

Given a perfect fluid (U,p = r ^-X r ,p = X r ) (or equivalently a local 
r-harmonic morphism) on a 4-dimensional spacetime (M, g) satisfying 
(fTUj) . we aim to prove that at some point where A ^ and for r^Owe 
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have either [/(In A) = (no expansion) or Q(X,Y) = (no rotation). 
The starting point is the observation that if the conjecture is true, 
then X(lnA), K(lnA) and Z(\nX) have to be basic (possibly all equal 
to zero). Moreover, a converse result holds (Proposition \5§ providing us 
with an equivalent form of the conjecture which turns out to be more 
tractable. Indeed, Equations f l26|) and fl27|) allow us to express X(lnA) 
and Y(lnX) as rational functions (with basic coefficients) of Z(lnX), 
provided that D = b 2 Z(\nX) 2 - b(c x - b 2 )Z(\nX) + (b x c 2 - b 2 a) is not 
zero, and then to produce a polynomial equation with basic coefficients 
in one variable, Z(\n A), which will be constrained to be basic (together 
with X(lnA) and Y{\nX)). So the proof will split up into two major 
cases (D / and D — 0) to be treated with independent methods. 
Nevertheless they have in common the following basic tools that will 
help us to conclude in each case. 

Lemma 2. If a function f on M satisfies a polynomial equation a n f n + 
... + ai/ + «o = where cti are all basic functions, then either f is basic 
or at = for all i. 

Proof. Iterate the derivative along V of the polynomial equation. At 
each step we have either V(f) = or an equation of smaller degree is 
satisfied. If / is not basic, i.e. V(f) ^ 0, then after n derivations we 
obtain ot n = 0. □ 

Lemma 3. If A ^ is a solution of the equation 

n 

where rji ? s are basic functions and r] io ^ 0, then either K(ln A) = or 
it exists j 7^ i such that p io = pj . 

Proof. Suppose V(ln A) ^ 0. Take (n — 1) times the ^/-derivative of 
the given equation and simplify it by V(lnA) and A Pj for every j ^ i . 
The resulted equation is VioY[j^i (Pio ~ Pj) = an d the conclusion 
follows. □ 

Let TZq = {2,3,4,y}. Recall that the conjecture is true for r 6 IZ Q . 

Proposition 5. Let (U,p = L ^-X r ,p = X r ) be a shear-free perfect fluid 
coupled with gravity on (M, g) and let {X, Y, Z} be horizontal vectors 
from a preferred orthonormal frame at a point of M. 

(i) If Z (In X) is basic and Z(lnX) ^ 0, then either the expansion or 
the rotation of the fluid vanishes. 

(ii) If X (In A) and Y(\nX) are basic and X(lnA) 2 + F(lnA) 2 ^ 0, 
then either the expansion or the rotation of the fluid vanishes. 

Proof. Since r 6 Hq are already settled, we may assume that r ^ IZ . 
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(i) By hypothesis Z(lnA) is basic, so Z(V(ln\)) = due to (ITT1) . 
Equation (TT8"j) becomes: 

(29) = ^^/3(Z) + (3-r)Z(lnA)y(lnA). 

Since Z(lnA) 7^ (by hypothesis), we have 

A 2(4-r) Q( Z \ 

Differentiating this equation along the vector X, implies: 

X{VQnX)) = ^((8 - 2r)/X(lnA) + X(f)) 

where the function / = — ( 3 _f^ lnA ) is basic. Inserting X(V(\ia. A)) from 
(ITS]) in the above equality gives us an equation with basic coefficients 
in X(lnA) and F(lnA): 

(31) (5 - r)fX(\n A) - (10 - 3r)fi(X, Y)YQn A) = /3(X) - X(/). 
Analogously we obtain 

(32) (10 - 3r)Q(X, Y)X(ln A) + (5 - r)/Y(ln A) = f3(Y) - Y(f). 

Since 0(Z) = — X 2 Q(X, Y)g([X, Y], Z) the discriminant of the linear 
system formed by (13~T1) and f|32|) is 

2 A 4 ^([x,y],z) 2 

j (3-r) 2 Z(lnA) 2 

Let us suppose that Q(X, Y) 7^ 0. Then A ^ (since r 7^ y ) and 
X(lnA) and F(lnA) are basic functions as solutions of a linear system 
with basic coefficients. 

By using (TBTJj) . the trace constraint (ITTj) becomes 

8-2,/ r-5 /3(Z) 2 3r-5 n . vv . 2 \ , 2r-5 u „ 2 



A = tt(X,Y) 2 (10-3r) 2 + (5 



A " ' 1,32(3 - r) 2 Z(ln A) 2 " ~^4~ fi(X ' F) J + "o" A 

- (2r ~ 5 j (r ~ 3) (X(ln A) 2 + y (In A) 2 + Z(ln A) 2 ) - ^ Scal^ = 0, 
12 

that is a linear equation in A 8_2r and A r ~ 2 with basic coefficients. Since 
the first two coefficients cannot cancel simultaneously, by applying 
Lemma [3] we conclude that [/(In A) = 0. 

(ii) Supposing X(lnA) 7^ 0, from ffTgj) we have 

A 2(4-r) p( X j + ( 10 _ 3r )y ( ln \)Q(x, Y) 



(33) K(lnA) 



4 (3-r)X(lnA) 

/3(X)+(10-3r)y(lnA)n(X,y) 



The function / = — 1 J (3- r )x(inA) ' * s basic according to our 
hypothesis. Suppose / 7^ (otherwise V(\n A) = and the proof ends). 
Analogously to the previous case (i), by differentiating ( 133]) along the 
vector Z and inserting Z(V^(lnA)) from (ITSi) we obtain 

(5-r)/Z(lnA)=|8(Z)-Z(/). 
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If r ^ 5, then Z(lnA) is basic. Inserting [/(In A) = ^-j—f in the trace 
constraint (fTTj) . we obtain again an equation of the form A 8_2r ■ basic + 
^F^A r ~ 2 + basic = 0; then, like in the case (i), we obtain [/(In A) = 0. 

If r = 5, let us suppose that Z(lnA) is not basic and show that 
this leads to a contradiction. From ( |2T|) we deduce that g(V xY, Z) = 
0, which inserted in (J25]) together with (133]) gives us: — /Z(lnA) 2 + 
^fQ(X,Y) 2 + 0(Z)Z(\nX) + basic term = 0. Derive along V this 
equation (and reinsert it into the result) to obtain 
/ 2 Z(lnA) 2 - /3(Z)fZ(lnX) + basic term = 0. According to Lemma | 
this implies / = 0, contradiction. □ 



5. The case D ^ 

Recall that D = b 2 Z(\nX) 2 - b(c x - b 2 )Z(\nX) + (b^ - 6 2 ci). In 
this case the key observation is that X(lnA), y(lnA) are rational func- 
tions with basic coefficients of Z(lnX) ^ 0. This allows, upon prop- 
agation of (126]) and (|27|) . to obtain a polynomial equation in Z(lnA) 
that leads us via Lemma [2] to the conclusion that Z(lnX) is basic and 
the conjecture is true according to Proposition [5j We mention that 
the proof makes use of the previously known fact [38J that the conjec- 
ture holds in the special case of aligned vorticity and acceleration, i.e. 
X(lnA) = F(lnA) = (see [SU] for a covariant proof). 

Assume Q ^ 0, D ^ and r £ TZ . From and (T27D we obtain 

(34) 

X(lnA) = i (bc 3 Z(\nX) 2 + (b 2 c 3 - b 3 c 2 + 6c )Z(lnA) + b 2 c - b c 2 ) , 
F(lnA) = J (-bb 3 Z(\nX) 2 + (6 3 ci - hc 3 - bb )Z(\nX) + b a - hc ) . 

Notice that we can also assume Z(lnX) ^ 0; otherwise from (13~4"]) we 
deduce that X(lnA), y(lnA) are both basic and either Proposition |5] 
applies (if one of them is not zero) or gradlnA G V (if they are both 
zero) and therefore U is irrotational (if A ^ constant), contradicting 
the assumption Q ^ 0, or [/(In A) = (if A = constant). 

By differentiating Equation (126]) along V, inserting second order 
derivatives from ( I18p and simplifying the result by means of (126 p and 
( ]2"7j) . we obtain 

(3 - r)V(ln A) {-&iX(ln A) - b 2 Y(\n A) - b 3 Z(ln A) - 2b } 
= ^{(10 - Sr)n(X,Y)( Cl + b 2 )X(\nX) 
(35) + [(10 - 3r)Q{X, Y)(c 2 - h) - bf3(Z)]Y(\n A) 

+ [(10 - 3r)Q(X, Y)c 3 - bp(X)]ZQn A) 
+ (10 - 3r)Q(X, Y)c - hf3(X) - b 2 /3(Y) - b 3 (3(Z)}, 
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while differentiating (127)) along V gives us 

(3 - r)V(ln A) {- Cl X(ln A) - c 2 F(ln A) - c 3 Z(\n A) - 2c } 

= ^{[(10 - 3r)Q(X,Y)(c 2 - 6 a ) + 6/3(Z)]X(ln A) 

(36) -(10-3r)fi(X,F)(c 1 + 6 2 )F(lnA) 

- [(10 - 3r)Q(X, Y)b 3 - b(3(X)}Z(\n A) 

- (10 - 3r)fi(X, Y)b - aP(X) - c 2 (3(Y) - c 3 /3(Z)}. 

Eliminating V(ln A) and A 8_2r from Equations (1331) and (13"6")) and in- 
serting X(lnA) and Y(lnX) from (131)) gives us a 6 th degree polynomial 
equation with basic coefficients in Z(lnA): V[Z(\nX)) = (see Ap- 
pendix for the explicit form). According to Lemma [U either Z(lnA) 
is basic or the coefficients of V are all vanishing. In the former case 
Proposition [S] applies and the conjecture is true, while in the latter we 
shall obtain a contradiction (except for three values of r). 

Let us suppose that Z(lnA) is not basic and that the coefficients of 
V are all zero, in particular the leading one: 

(37) b\ + cj - £5§ (b 3 /3(X) + c 3 /3(Y)) = 0. 

Lemma 4. and Z(lnA) is not basic, then b 2 + C\ = 

and bi — c 2 = 0. 

Proof. Substitute X(ln A), K(lnA) from (J31} and X(Z(\n A)), Y(Z(\nX)) 
from (1211) . in Equation (I22p to obtain a 7 th degree polynomial equation 
with basic coefficients in Z(lnA) whose leading coefficient b 6 (b 2 + c±). 
Since Z(lnA) is not basic, b 2 + c\ must vanish according to Lemma [2j 
Analogously, substituting X(lnA), F(lnA) from ( 134]) in Equation 
(12 ip. then summing with the corresponding relation with X and Y 
permuted, we obtain again a 7 th degree polynomial equation in Z(ln A) 
with basic coefficients, the leading one being b 6 (bi — c 2 ). □ 

This Lemma allows us to eliminate bi and b 2 in the following com- 
putations. According to ( 131)) and (133)) . V(lnA) = ^-j— <3, where Q = 
(b 2 c 3 Z(\nX) 3 + ...)/((r - 3)b 2 b 3 Z(\nX) 3 + ...) is a rational function in 
Z(lnA) with basic coefficients. Plugging this into ( 12 5 j) together with 
X(lnA) and Y(lnX) from (I34p gives us: 

A 8_2r b 6 u(r)c 3 Z (In A) 9 + . . . (lower degree terms) 

4 6 6 t> (r)Q(X, Y) 2 c 3 Z(ln A) 7 + . . . (lower degree terms) ' 
where c 3 := c 3 - ^(F), u(r) = (^(^-s) and w(r) = ^-y . 

Suppose r 7^ 5/2 to have u(r) 7^ (notice that t>(r) 7^ 0). Replacing 
V(lnA) and A§ 4 r with their (rational) expression in terms of Z(lnA) 
in the identity V(^-) = (8 - 2r)^-V(\nX) and using (JTHD gives 
a polynomial equation with basic coefficients in Z(ln A) whose leading 
term vanishes if and only if 63 = or c 3 = 0. A similar argument, 
starting from ( 136)) . imposes c 3 = or b 3 = 0, where b 3 := b 3 — ^|/3(X). 
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Due to (157)) we have to consider only the following two sub-cases. 

5.1. Sub-case 6 3 = C3 = 0. Returning to the polynomial equation 
V = mentioned above, we obtain b (3(X) + c Q (3(Y) = from the 5 th 
degree coefficient term. Inserting it in the 4 th and 3 rd degree coefficient 
terms we obtain respectively 

2{bl + cg)(10 - 3r)Q(X, Y) + 3c 2 (b f3(Y) - c f3(X)) = 0, 

ci [7{b 2 + c 2 )(10 - 3r)n(X, Y) + 10c 2 (b f3(Y) - c (3(X))} = 0. 

If ci 7^ 0, the system (J3SJ) implies b = c = 0; therefore X(lnA) = 
K(lnA) = 0, that is we are in the known case of aligned vorticity 
and acceleration where the conjecture is true |38j, so V(lnA) = that 
contradicts our assumption that Z(lnA) is not basic. 

If c\ = 0, the coefficient of the 2 rd degree term gives us 



^Ef (K + c 2 )(W - 3r)n(X, Y) + 4c 2 (b p(Y) - c /3(X)) 



0. 



where we have used f3(Z) = — ^-4) ° 2 ' consec L uence of the fact es- 
tablished in Lemma HI So either c 2 = (and, from (1381) we obtain 
again b = c = 0, so X(lnA) = F(lnA) = 0) or c 2 7^ and the system 
formed by the above equation and the first Equation (155)) implies again 
bo = Cq = (so X(\n\) = F(lnA) = 0) except for r 
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5.2. Sub-case 63 = C3 = 0. In this case the system formed by the 
cancellation of coefficients of V can be solved in c± and c 2 and, after a 
case by case analysis, we obtain (except for r = -r) 

X(lnA) = ( 10 _3r)n(x,y) > y(lnA) = — (io-fr)n|x,y) • 

Plugging this into (1T8)) we see that either X(ln A), K (In A) are both zero 
(the case of aligned vorticity and acceleration) or V(lnA) = 0, which 
contradicts the assumption that Z(lnA) is not basic. 

In conclusion, if D 7^ 0, the conjecture is true for r ^ {y, |, y}. 

6. The case D = 

In this case we deal essentially with the situation when the vorticity 
is orthogonal to the acceleration, Z(hxX) = 0. Since the compatibility 
conditions ( 126)) and ( )27j) become trivial, the remaining one, ( 128)) . comes 
into play. By successive propagation of this condition we can eliminate 
all terms in X(lnA) and F(lnA) and obtain a quadratic equation in 
V(V(lnA)) and V(\n\) 2 with coefficients involving only A and basic 
quantities. This quadratic equation conserves its form upon iterated 
^-derivations, so the terms in V^(lnA) can be eliminated and finally it 
results an equation for which Lemma |5J applies to conclude that the 
fluid must be expansion free (once the rotation is non-zero). 

Assume that D = and r ^ IZq. According to Lemma EJ either 6 = 
(that is Q = and the proof ends) or Z(ln A) is basic. If Z(\n A) is basic 
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and Z(lnA) 7^ 0, Proposition applies and we obtain the conclusion. 
The remaining case is Z(lnA) = 0. 

For the rest of this section, let us suppose moreover that Z(ln A) = 
and ft(X, Y) ^ 0. 

From (US) we see that 0(Z) = that is g([X, Y], Z) = (recall that 
f3(Z) = -X 2 g{[X,Y],Z)Q{X,Y)). Equation ([26]) reduces to 

(39) X 2 g(V x Y,Z)X(\nX) + X 2 g([Z,Y],Y)Y(\nX) = ggzafeyy , 

while Equation ( ]2"7|) reduces to 



(40) A 2 s([Z,X],X)X(lnA) + A 2 <7(VxY, Z)y(ln A) = - (10 _ 3r) o(x,Y) ■ 

Notice that the system formed by (139]) and (1401 is dependent (L> = 0). 

Therefore two situations have to be considered: in ( 13 9 p at least one 
of the (basic) coefficients of X(lnA) and F(lnA) is not vanishing, or 
both are zero. 

6.1. Sub-case g(V x Y, Z) 2 + g([Z,Y},Y) 2 ^ 0. Suppose g(V x Y,Z) ^ 
(the other choice is similar). According to (139]) . X(lnA) is a linear 
function of F(lnA) with basic coefficients: X(lnA) = aF(lnA) + b. 
Taking the ^-derivative of this relation and employing ( II 8p we obtain: 

VQn\) = i^{P(aX-Y) + (W-3rMX,Y)[(a 2 + l)Y(\nX) + ab]}, 

where we supposed b ^ (otherwise the conclusion is immediate). 
Plugging this into Equation (1231) we obtain a third degree polynomial in 
y(lnA) with basic coefficients. According to Lemma [2] either y(lnA) 
is basic or the coefficients are all vanishing. If y(lnA) is basic then 
X(lnA) is basic and we can apply Proposition [5] to conclude. If this is 
not the case, then the leading coefficient (a 4 — 1) (10 — 3r) (3 — r) (5 — r — 
A(r))Q(X, Y) must vanish, so r = 5/2 ( not covered by our proof) 
or a = ±1. But in the last case, the coefficient second degree term is 



8(r-4) 
r-2 



(10 — 3r)(3 — r)bQ(X, Y) and it cannot be zero, contradiction. 



6.2. Sub-case g(V x Y, Z) = and g([Z,Y],Y) = 0. Inserting this in 
( 140]) we obtain that either X(lnA) is basic, or g([Z, X], X) = 0. 

Firstly, let us suppose that g([Z, X],X) 7^ (so X(lnA) is basic). 
From ( ITS]) we have: 
(41) 

^f^{(5{X) + (10 - 3r)F (In A)ft(X, Y)} + (3 - r)X(ln X)V(ln A) = 0. 

If X(lnA) = 0, then F(lnA) = — (lo-^npc y) ^ s basic, so either Propo- 
sition ED applies (if F(lnA) 7^ 0), or gradlnA e V (if F(lnA) = 0). In 
the latter case U is hypersurface orthogonal (so irrotational, a contra- 
diction) if A 7^ constant, or [/(In A) = if A is constant. 
If X(lnA) 7^ 0, then eliminating V(lnA) and from (J23D and (jUJ) 
we obtain a 3 rd degree polynomial equation in Y (In A) with basic coeffi- 
cients whose leading coefficient is (10 — 3r)(3 — r)(5 — r — A(r))ft(X, Y). 
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Again Lemma [2] imposes that either K(lnA) is basic (so Proposition 
applies) or the leading coefficient is zero, i.e. r = 5/2 (excepted case). 

Secondly, let us suppose that g([Z,X],X) = 0. This remaining sub- 
case can be resumed by the following assumptions: 

• Z(lnA) = 0; 

• g(V x Y,Z)=g(V Y X,Z) = 0; 

• g([Z,X],X) = g([Z,Y],Y) = (=► Z(Q(X,Y)) = 0, divZ = 0). 

In this last sub-case Z itself is shear-free and divergence-free Q. More- 
over its 3-dimensional orthogonal distribution is integrable (so Z is 
aligned with a gradient vector) and minimal. It also follows from the 
above conditions that Hessi n A(Z, X) = Hessi n A(Z, Y) = 0, (f* Kic N (Z, X) = 
cp* Ric N (Z, Y) = and d(3(Z, X) = d(3(Z, Y) = 0. 

Since, in this case, Hessi nA (Z, Z) = -X(ln A) 2 -F(ln A) 2 +5iX(ln A) + 
5 2 Y(\nX) + A~ 2 £/(lnA) 2 , Equation (TJOD becomes: 

(1 -£(r))(X(lnA) 2 + y(hiA) 2 ) -SiX(ln\) -5 2 F(lnA) 

(42) - iA- 2 f/(ln A) 2 + |A^ 2 - »ffi\*-»n(X, Yf 
- ^ Ric N (Z, Z) + ^ Scal^) = 0, 

where we introduced the notations 8\ = \ 2 g([Z, X], Z) and 5 2 = \ 2 g([Z, Y],Z). 

Recall that (5{X) = Y{Q(X, Y))+5 2 Q{X, Y) and (3{Y) = -X{Q{X, Y)) — 
5 l Q{X,Y) (see the proof of Lemma [T]). Divide (J2HD by &{X,Y), mul- 
tiply (j4"2|) with (10 — 3r) and take their sum to obtain: 

(C(r) + (10 - 3r)(l - B(r)) [X(lnA) 2 + F(lnA) 2 ] 

( 43 ) + 30_7r A - 2[/(ln A) 2 + 11^30 y~2 + 60+44r-21r" A 8-2r fl(X> y) 2 

+ Scal^ -»?V Ric^(Z, Z) + = 0. 

From ([18]) we deduce the following propagation equations: 

V (X(ln A) 2 + y(ln A) 2 ) = 2{^ [/3(X)X(lnA) + /3(F)r(lnA)] 
+ (3-r)V(InA) (X(ln A) 2 + YQn A) 2 ) }, 

V (/3(F)X(ln A) - /3(X)F (In A)) 

= (10 - 3r)Q{X, Y)^ [P(X)XQn A) + ${Y)Y (In A)] 
+ (3 - r)V(mA) (/3(F)X(lnA) - /3(X)F(lnA)) , 



We are close to the situation in which there exists a Killing vector parallel to the 
vorticity, a notorious case for difficulty despite the drastic simplifications it brings 
about, see [TP] . 
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V (P(X)X(]n\) + /3(Y)Y(\nX)) 

= ^ {/3(A) 2 + /3(F) 2 - (10 - 3r)n(X, Y) [/3(Y)X(ln A) - 0(X)Y(\n A)] } 
+ (3 - r)V(ln A) (/3(A)X(ln A) + (In A)) . 

The "cyclic" property of the ^/-derivatives of the linear term P(Y)X(ln A) — 
f3(X)Y (In A) allows us to eliminate it by taking the derivative of Equa- 
tion ( 143|) twice along V. 

The first ^/-derivative of Equation ( T43]) along V (simplified by rein- 
serting p(Y)X(\n\) - p(X)Y(\n\) from (03]) in the result) reads: 

(44) 

^"Vr' 30 ^ [0(*)*O* A) + (In A)] 

+ V (In A) { ^-3)V-_4)(7r-10) [x(ln A)2 + y (ln A)2 ] 
2-r(2r-5) x r-2 , r(3r-14)(7r-20) x 8~2rry y y\2 

H 3 — a h — A ^tA,r; 

- (r - 3) (|S| Seal" +2!fc*V Hfc"(Z, Z) - } = 0. 

Now, taking the derivative along V of Equation (j4"4"|) . then replacing 
/3(X)X(lnA)+/3(F)y(lnA) from (|SD, /3(y)X(ln A)-/3(A)F(ln A) from 
(T4"3"j) and also X(lnA) 2 + F(lnA) 2 from the trace constraint (TTTj) . gives 
us an equation of the following form: 

(45) 

a(r)V(VQn A)) 2 + 0(r)V(V(ln A))V(ln A) 2 + 7 (r)y(ln A) 4 

+ \/(lnA) 2 P!(A 8 - 2r ,A r - 2 ) + \/(\/(lnA))P 2 (A 8 ~ 2r ,A r - 2 ) + P 3 (A 8 ~ 2r ,A r ~ 2 ) = 0, 

where Pj are polynomial expressions with basic coefficients (see Appen- 
dix for the explicit form) and we supposed r ^ {§' fi> t) ^° nave wen 
defined expressions in r. 

From ( I17j) we obtain the 3 rd order ^/-derivative of In A in the form: 
(46) 

V(V(V(\nX))) = V(lnX) ( gi (r)V(V(\n A)) + q 2 (r)V(\n A) 2 + Po(A 8 " 2r , A r - 2 )) , 

where we have used Equation (1441) to eliminate f3 (X) X (In X)+ (3 (Y)Y (In X) 
and again f fT7|) to eliminate X(lnA) 2 + F(ln A) 2 (see Appendix for the 
explicit form). This constraint on V(V(V(ln A))) implies that Equa- 
tion (|45p conserves its general form when it is derived along V (after 
simplification of V(ln A), if not zero). More precisely, the propagation 
of (I4"5"|) . after inserting (T4"6"l) and dividing by V(\nX), is 

(47) 

a {1) (r)V(V(\n A)) 2 + (3 m (r)V(V(\n X))V(\n A) 2 + 7 (1) (r)K(ln A) 4 

+ y(lnA) 2 Pf ) (A 8 - 2r ,A r - 2 ) + y(\/(lnA))P 2 (1) (A 8 - 2r ,A r ' 2 ) + P 3 (1) (A 8 - 2r ,A r " 2 ) = 0, 
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with 

a W( r ) = a(r) qi (r) + 2(3{r), P x (1) = P[ + q 2 {r)P 2 - q 1 (r)P 1 + 0(r)P o , 
/3 (1) (r) = 2a(r)q 2 (r) + 4 7 (r), P 2 (1) = P' 2 + 2 Pi + 2a(r)P , 
7 W(r) = P(r)q 2 (r) - *y(r) qi (r), P 3 {1) = Pi - qi (r)P 3 + P P 2 , 
where P[ = AfJ. 

Iterating the V^-derivative of ( 14 5 p shows that a^ k \ (3^, 7^ form geo- 
metric progressions with the same common ratio (or (a^,^ 1 ^' 1 ') is 
an eigenvector of the recurrence matrix). Therefore, from the first four 
^-derivatives of (145]) we obtain an overdetermined linear system in 
^(lnA) 2 and V(V(ln\)). Its compatibility condition (zero determi- 
nant) leads us to a 4 th degree polynomial equation P(A 8_2r , A r ~ 2 ) = 
having basic coefficients: 

where r/i and r\ 2 are rational functions with common zeros r = 0, r = — 
and r = y (these are common zeros for all 4 th degree coefficient terms) . 

Recall that r 7^ by the conjecture hypothesis and that r / | as 
we have already assumed. Suppose moreover r ^ |. If rji(r) = 0, we 
employ Lemma [3] with reference coefficient 772. Since the corresponding 
exponent of A is different from the others exponents in V, we conclude 
that V(lnA) = 0. If rjx (r) ^ we apply Lemma |3] with reference 
coefficient T]i to deduce that either V(lnA) = or r £ {y, y,6}. But 
in the latter case we check that Lemma [3] with reference coefficient 772 
leads us to V(ln A) = 0. 

In conclusion, if D = 0, the conjecture is true for r ^ {|, fj, y, -j}. 

7. Conclusion 

We have proved the following 

Theorem 1. Consider a shear-free perfect fluid solution of Einstein's 
field equations where the fluid pressure satisfies a barotropic equation 
of state of the form p = wp with 

w G R \ {-§, -i, -i, |} . 

Then the fluid is either non-rotating or non- expanding. 

The conclusion of the theorem is in fact sharper. On one hand, in 
the expanding and non-rotating case, the acceleration vanishes (and 
grad^ In A = 0) and the spacetime must be in this warped prod- 

uct of an interval with a constant curvature 3-manifold, that is a FRW 
model. This was proved in [9] by analysing the complete list of space- 
times that admit a shear-free non-rotating fluid, classified in terms of 
the Weyl curvature symmetries. Another, direct proof can be done by 
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starting with compatibility equations as ( I2"3"j) . (I2"5"j) and the propaga- 
tion of (T2TT) . and showing that grad^lnA = (this requires a separate 
analysis for r £ {3,4, |}). 

On the other hand, in the non-expanding and rotating case, the 
spacetime must be stationary since the fundamental vector V is Killing. 

Thus, adopting the dual perspective of r-harmonic morphisms on 
spacetimes satisfying ([10]) . the above result states that they are either 
of warped product type or of Killing type (with six exceptions for r). 

Concerning the exceptional values of w, we notice that the two posi- 
tive ones are in the physical regime defined by a speed of sound between 
and 1, and that all excepted fluids are satisfying the strong energy 
condition p + 3p > 0. We expect that the proof can be adapted for the 
remaining values of w as for the already known cases w £ {0,±|,|} 
that also appeared to be exceptional in our proof. To what extent it 
can also be adapted for non-linear equations of state is the next step 
to consider. We postpone these tasks for a future work. 

8. Appendix 

8.1. Notations and conventions. Throughout the paper (M,g) de- 
notes a connected time-oriented Lorentz 4-dimensional manifold with 
metric signature (—,+,+, +). We denote by £ the Lie derivative, 
by V the Levi-Civita connection of (M,g), and we use the following 
sign conventions for the curvature tensor field R(X, Y)Z = V ^VyZ — 
VyVx-Z — V[x,y]Z, and Af = div(grad/) for the Laplacian on func- 
tions. Hessj = Vdf is the Hessian, while Ric and Seal denote the Ricci 
and scalar curvature of (M,g), respectively. 

Given a 2-form Q on (M,g), the interior product with a vector X 
is = Q{X, •) and its co-differential is defined in analogy with Rie- 
mannian case by 5Q(X) = — J2t=i £i(V ei f2)(ej, X), where {ei} i= j^ is an 
orthonormal frame with g[e^ e^) = £j = ±1. The metric g on M induces 
a (pointwise) metric on the bundle of p-covariant tensors (or p- forms) on 

M, defined by: (A, B) = ^ E? lr .. )ip =i ■ ■ ■ • • • > e i P ) B ( e h > • • • > e %,)> 

that provides us with the norm |^4| of such object. The notation if* A 
refers to the usual pullback of A by the mapping ip. 

Along the paper we extensively used Mathematica ([39]) to simplify 
the coefficients and to solve polynomial equations. Many of the formu- 
lae involved in the final argument in both principal cases are very long 
and were not included in the main text. We reproduce here some of 
them. Full version is available in separate files. 

8.2. Details for the case D ^ 0. When h\ = ci and b 2 = —C\ the 
polynomial we employed is V(Z(ln A)) = Y^=o Coeffj Z(lnA)*, where 

Coeff 6 = -6 4 [b(hP(X) + c 3 (3{Y)) - (6| + c 2 3 )n(X, F)(10 - 3r)] , 
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Coeff 5 = -b 3 [2b 2 (b f3(X) + c /3(Y)) - 2b((2b 3Cl - c 2 c 3 )f3(X) + (b 3 c 2 + 
2 Cl c 3 )P(Y)) + 3 ( Cl (&§ + 4) - b(b b 3 + c c 3 )) (10 - 3r)Q(X, Y)] , 

Coeff 4 = -b 2 [ - (10 - 3r) (3c? + c|) (b 2 + c 2 )Q{X, Y) + b 2 (b 3 (3{Z)c - 
& (9/3(X)c 1 + 3/3(F)c 2 + /3(Z)c 3 )-26^(X,r)(10-3r) + co(-9/3(F)c 1 + 
3/3(X)c 2 -2(10-3r)con(X,r)))-6(62/3(Z)c 2 -26 3 c 1 (3/3(X)c 1 +3/3(r)c 2 + 
56 n(X, y) (10 - 3r)) + c 3 ( - 6/3(y )c? + 6/3(X)cic 2 +/3(Z)c 2 c 3 - 10(10 - 
3r) CoCl fi(X,y)))], 

Coeff 3 = b[- Cl (c 2 + c 2 ) (b 2 + c 2 )Q{X, y)(10 - 3r) - 2&(& 2 /3(Z) Cl c 2 - 
b 3 (2f3(X)c\+ (3c?+c 2 ) (f3(Y)c 2 +2b Q(X, Y ) (10-3r))) +c 3 (-2/5(y)c?+ 
f3(X)(3cjc 2 + 4) + /3(Z)c lC2 c 3 - Wco4tt(X,Y) - 20c c 2 ft(X, y) + 
6c (3c? + 4)Q(X,Y)r)) + b 2 (2b 3 p(Z)(2c oCl + b c 2 ) - 26 (8/3(X)c? + 
5/3(y)cic 2 +/3(X)c2+2/3(Z)cic 3 )-76gcin(X,y)(10-3r)+co(lO/3(X)c lC2 - 
2/3(y) (8c? + c|) + 2f3(Z)c 2 c 3 - 7c oCl fi(X, y)(10 - 3r)))] , 

Coeff 2 = 6[6 2 /3(Z) (6g + c 2 )c 2 + 6( - 26 3 /3(Z) (26 oCl c 2 + c (3c? + 4)) + 
26 (/3(X)(7c? + 3c lC 2 ) + (3c? + c 2 2 )(2f3(Y)c 2 + P(Z)c 3 )) + 36 2 (3c? + 
c|)^(X,y)(10 - 3r) + c (2/3(y)(7c? + 3c lC |) - 4c 2 (/3(X) (3c? + c 2 ) + 
/3(Z) Cl c 3 ) +3c (3c? + C 2)(10-3r)fi(X,y))) + (c? + 4) \b 2 c 2 f3(Z) + 
6 3 (P(X) (-4+4) -2ci (/3(y )c 2 +36 fi(X, y ) (10-3r))) +c 3 (2/3(X) Cl c 2 + 
/3(y) (-4 + 4) + P(Z)c 2 c 3 - 6c oCl Q(X, 30(10 - 3r)))] , 

Coeff! = -26 2 (6g + c 2 ) Cl c 2 /3(Z) + (c? + c 2 ) 2 (& & 3 + c c 3 )fi(X, y)(10 - 
3r) + 6(c? + 4) (2b 3 P{Z){2c oCl + b c 2 ) - 26 ci(3/3(X)ci + 3f3(Y)c 2 + 
2/3(Z)c 3 )-5^ Cl fi(X,y)(10-3r)+c (-6c?/3(y)+6cic 2 /3(X)+2c 2 c 3 /3(Z)- 
5c cifi(X,y)(10-3r))), 

Coef? Q = (c 2 + c 2 2 )[-b 3 c (c 2 1 + 4)f3(Z)+b (c 2 + 4)((3(X)c 1 + (3(Y)c 2 + 
0(Z)c 3 ) + c {bf3(Z)c c 2 + (c? + 4) ((3(Y) Cl - (3(X)c 2 + c Q(X, Y)(10 - 
3r))) + b 2 (bp(Z)c 2 + (c? + 4)Q(X, y)(10 - 3r))] . 



8.3. Details for the case D = 0. In Equation (Hoi) . 

7(r) = a(r) fflz3gg^a ; 

and 

Pi(A 8 ~ 2r , \ T - 2 ) =h\ 8 - 2r + k 2 X 6 ~ r + A; 3 A 16 " 4r , 
P 2 (A 8 - 2r , A r - 2 ) =£ 1 X 8 - 2r + £ 2 X 6 ~ r + £ 3 A 16 " 4r , 

P 3 (A 8 " 2r , A r ~ 2 ) = (r " 3)( 2 T (r 4 S ) lr " 3 ° ) ^(X y) 2 A 16 " 4r (rmA 8 -^ + m 2 A 6 ~ r + m 3 A 16 " 4r ) 
where 

7 _ ( 3 -r) 2 (7200-5940r+1760r 2 -257r 3 +21r 4 ) q ijy (r-3)(r-5)(19r-54) a 
Kl ~~ 3(r-2)(2r-5) (llr-30) SCal llr-30 A ^> 



-21600+39840r-26944r 2 +8588r 3 -1312r 4 +78r 5 
3(r-2) (llr-30) 
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U 12960000-49248000r+70419600r 2 -53293200r 3 +24116984r 4 -6806516r 5 +1184369r 6 -117200r 7 +5082r 8 
3 24(r-2) 2 (2r-5)(4r-15)(llr-30) 

t^- ^X-f Scal"+(3-r)A„ 

p _ 2(-120+164r-68r 2 +9r 3 ) 
12 ~ 3(r-2) ' 

p _ 144000-422400r+439960r 2 -224068r 3 +60574T- 4 -8353r 5 +462r 6 p/ V V\ 2 
*-3 — 24(r-2) 2 (2r-5)(4r-15) I ) , 

_ 3(10-3r) 2 q ]N , 10-3r a i /3(X) 2 +/3(y) 2 

'"1 — 2(2r-5)(4r-15) OL,d ^4r-15 iiz ^ Q(X,y) 2 ' 

_ _ r(r-6)(10-3r) 2 
m 2 ~~ 6(r-2)(r-3)(4r-15) ' 

_ (10-3r) 2 (30-23r+5r 2 ) ry y y\2 

" 4 3 — 4(r~3)(r-2)(2r-5)(4r-15)"v^-5 2 J 5 

where A z := 2±^fV Ric"(Z, Z) - ^gg^. 
In Equation (J15I) . 

_ /„\ 91r 2 -547r+810 „ /„\ (3r-ll)(29r 2 -168r+240) 

and 

p (\8-2r \ '/" 2 \ _ / ( 3 -r) 2 (r-2)(10-3r) q ,7V (3-r)(r-2)(2r-5) a \ )8-2r 

P)l A J - I (r-4)(llr-30) &Cal + 2(r-4)(llr-30) A ^ J A 

A 6 



(2r-5)(10-3r)(17r 2 -92r+120) ,6-r 



6(r-4)(llr-30) 

, r(10-3r)(7r-20) x 16-4rQ/ y V\ 2 
~^ 8(r-4)(llr-30) A "^i 1 i 

or, in short notation, P (A 8 ~ 2r , A r ~ 2 ) = niA 8 ~ 2r + n 2 A 6 ~ r + n 3 A 16 " 4r . 

At the level of coefficients, the recurrence relations involved in the 
iterated K-derivatives of Equation ( H7|) read: 

= -KifiaV-VipV) = -K{r)^- l \^ = -K^r)^- 1 ^, j > 2 

where K(r) := 2b(j ^-m^ anc ^ ^ or 3 — ^-1 

kf = (6 - r - qi(r))ki j ' 1] + ®(r)^ _1) + /3^- 1 )(r)n 2 , 

fc^ = (16 -At- qi{r))k { i~ 1] + g 2 (r)4 i_1) + /3 (i_1) (^)^3, 

4° = (6 - r)4 i_1) + 2k [ t l) + 2a^'- 1 )(r)n 2 , 

4 j) = (16 - 4r)4'" 1) + 24 J_1) + 2a^'- 1 )(r)n 3 , 

where fc^ = hi, if = £j, for all i, = a and (3^ = j3 (the analogous 
recurrence relations for ki and tf 1 have not been used). 

For P 3 (j) , the recurrence relation gets slightly more complicated since 

pW( A 8-2r j A r-2) = A 16-4r [ m (l) A 8-2r + m « A 6-r + m (D A 16-4r 

+ TO W A 2-4 + TO (l) A r-2 + m (l)] ) 

where we focus on the following coefficients: 
m!p = (22 - 5r — <7i(r))m 2 + £ 2 n 3 + 4^2, 
= (32 — 8r — gi(r))m3 + ^3^3, 



.CO 



with m % := (r - 3)( 24 ^ ) lr - 30) ^(X,F) 2 m t . 

Then, for j > 2 , P3 will have the same form and we have 
= (22 - 5r - gi(r))7n^ -1) + ^J'" 1 ^ + 4 i_1)n 2, 
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mi' = (32 - 8r - qi{r))w>i iJ + 4 J >n 3 , 
mf = (12 - 2r - ?i(r))mj' _1) + 

Taking linear combinations of the ^-derivatives of Equation (T4"T|) . we 
formed the overdetermined linear system in K(lnA) 2 and V^V^lnA)): 

(Pi j) + K{r)P[ J ^ ^(lnA) 2 +(P 2 (i) + #(r).P 2 (5 '~ 1) ) V(V(\n A))+P 3 (j) + 

K(r)P!f ^ = 0, where j = 2,3,4. The determinant of its coefficients 
matrix (that must vanish) is, after dividing by \^ 8 ~ 2r \ a polynomial 
P(A 8_2r , A r_2 ) with basic coefficients. For instance the coefficient of 

A 4(8-2r) ig 





1.(2) 
3 


+ iT(r)4 1} 


/ 2) 


+ K(r)l { 3 ] 


+ iiT(r] 


(i) 
m 3 


n(x,Y) 8 


1.(3) 
3 


+ K(r)kf ] 


*3 




+ K(r] 


(2) 

m 3 




1.(4) 
ft 3 


+ K(r)4 3) 


«(4) 
*3 


+ K(r)tf 


mj 4 ' + if(r' 


(3) 
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